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In this paper we extend the concept of separatrix reconnection into chaotic regimes. We show that 
even under chaotic conditions one can still understand abrupt jumps of diffusive-like processes in 
the relevant phase-space in terms of relatively smooth realignments of stable and unstable manifolds 
of unstable fixed points. 
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Two dimensional nonmonotonic conservative maps are 
recognized to be of relevance in modeling a number of 
nonlinear systems, as for instance laser acceleration of 
charged particles and the nonlinear flow of mag- 

netic field lines in fusion machines like tokamaks and 
others |t]||. As opposed to the more traditional mono- 
tonic versions, nonmonotonic maps are characterized by 
frequency curves that are not monotonic functions of the 
action variable. In laser accelerators, nonmonotonicity 
arises as a result of the relativistic mass variation of the 
accelerating particles j|; in tokamaks, it arises as a re- 
sult of the geometrical peculiarities of the relevant back- 
ground magnetic fields. In any case, nonmonotonicity has 
strong influence on the types of bifurcations that can oc- 
cur in the associated nonlinear dynamics. Indeed, mono- 
tonic maps typically allow only period doubling bifur- 
cations of fixed points but if one adds nonmonotonicity, 
tangent bifurcations involving pairs of elliptic and hyper- 
bolic fixed points also become possible in the appropriate 
phase-space. 

Period doubling cascades of periodic orbits generally 
precede a transition to chaotic regimes of these orbits, 
but tangent bifurcations bear no direct relationship to 
nonintegrability. Indeed, it has been argued that the 
most noticeable effects of a tangent bifurcation are to be 
seen while the system lies in predominantly integrable 
regimes. Although we will show that this sort of argu- 
ment is questionable, consider the process depicted in the 
integrable case of Fig. (|I|) Q . Two chains of fixed points 
undergo tangent bifurcation, or actually an inverse tan- 
gent bifurcation, starting from the leftmost panel down 
to the rightmost. Before the bifurcation itself where el- 
liptic fixed points collapse against hyperbolic points, the 
separatrices defining the upper chain undergo a recon- 
necting process with those defining the lower chain - this 
is seen in Fig.(|l|b). It is precisely due to the smooth- 
ness caused by integrability that the reconnection can be 
seen so clearly. This is why reconnection is thought to 
be of relevance only in integrable cases. In contrast, the 
process is generally regarded as of little significance in 
chaotic regimes because in those situations all separatri- 



ces - which shall be correctly called stable and unstable 
manifolds then - would be already interlaced with little 
global response as relevant control parameters are var- 
ied. Speaking in more precise terms, effects associated 
with reconnections are thought to be unobservable when 
the elliptic fixed points of the reconnecting chains un- 
dergo full cascades of period doublings, before any sort of 
mutual contact of the relevant manifolds takes place. 

While it is true that reconnections are not easily visu- 
alized in chaotic regimes, it is our purpose here to show 
that their effect can still be quite appreciable. What 
happens is that even in chaotic regimes the unstable or 
stable manifolds of originally hyperbolic points may still 
make a transition from a situation where their mutual 
crossings are absent - or, actually, relatively infrequent 
- to a situation where the mutual crossings become very 
frequent. We shall illustrate the process with a set of fig- 
ures from which one will be able to see that this change 
in topology is a reminiscent of the corresponding behav- 
ior of regular regimes. The macroscopic result of this 
type of transition is that as soon as the crossings become 
frequent, stochastic diffusion undergoes an abrupt jump 
from slower to faster rates. 

From this point on, the discussion relies on more tech- 
nical grounds. Therefore we introduce here the model 
map we shall be working with. The map is called the 
nonmonotonic twist map and reads: 

Pn+l =Pn - ksm(f) n , (j) n+ i =((>„ + f(p n +l), (1) 

where (p, </)) is a pair of discrete canonical variables - p 
representing an action and <f> a 27r-periodic angular coor- 
dinate - and where the map itself is totally simpletic given 
that the left-hand-side of the second Eq.(]l|) depends on 
Pn+i instead of p n . Function f(p) is of foremost impor- 
tance here. It is in fact a measure of the frequency with 
which the discrete orbits move on the (p, <fi) phase-space. 
In standard monotonic maps it reads f(p) — p, but as 
we wish to incorporate nonmonotonic features we add a 
quadratic term such that / becomes f(p) — p — ap 2 as in 
Rcf. jjj. In this case one has effectively a nonmonotonic 
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frequency curve with maximum located at p = l/(2a), 
where df(p)/dp — 0. The map (0) has several fami- 
lies of fixed points. Let us focus here on the first or- 
der family (period one orbits) which is characterized by 
Pn+i = Pn and <j) n +i — 4> n + 2mTr, with m as an arbitrary 
positive or negative integer; we shall refer to the fixed 
points as (Pmi^m)- Eq.(W) informs that the fixed points 
are located at §* m = 0, 7r andp* t — [1±-^1 — 8mira/(2a)]. 
In this case of period one orbits, let us have a brief look 
at the distribution of the various fixed points over the 
phase-space. For m = one has four points located at 
p* Q = 0, 1/a. If m > the fixed points lie in the fi- 
nite interval < p < 1/a, and if m < the points lie 
in any of the intervals — oo<p<0orl/a<p< +oo. 
The existence of points located in the finite interval must 
satisfy the condition m < l/(87ra), so it may well hap- 
pen that no fixed point can be actually found there if a 
is large enough. What happens is that as a grows from 
some small value, all the fixed points originally located in 
the sub-interval < p < l/(2a), collapse, at p = l/(2a), 
against the corresponding points originally located in the 
sub-interval l/(2a) < p < 1/a via a sequence of inverse 
tangent bifurcations. Meanwhile, the points placed ex- 
ternally to the interval (0, 1/a), as well as the points of 
the pair m = 0, simply approach each other but never 
touch. We point out that although this latter two types 
of fixed points never undergo inverse tangent bifurcation, 
their manifolds can naturally undergo the reconnection 
processes. In previous works |l[] the m = 1 case has been 
investigated. It has been shown that when chaos is ab- 
sent, in the sense that elliptic points of the various chains 
have not yet period doubled to chaos, initial conditions 
at negative values of p do not move up to positive values 
unless the separatrices of the two chains corresponding 
to the to = 1 resonances touch each other. In addition it 
has been argued that for those situations where the ellip- 
tic fixed points have already undergone full cascades of 
period doublings, no detectable difference in the global 
aspects of the dynamics should be observed as arising 
from a possible contact of separatrices. We now pro- 
ceed to show that this is not quite exact; even under 
chaotic conditions, some noticeable effects resulting from 
manifold reconnections can be in fact observed. Specifi- 
cally we shall show that after what could be best called 
a reconnection-like process, diffusion makes a somewhat 
abrupt jump from lower to higher values. 

Let us focus the discussion on the m = case because 
this resonance is the largest one in the system. The lin- 
ear stability of the fixed points can be examined from the 
characteristic equation A = [(2 — «) ± \/{k— 2) 2 — 4]/2, 
with k = cos(</>*)(l — 2ap*)k. A is the eigenvalue of the 
linearized map; if complex (purely real) the correspond- 
ing fixed point is unstable (stable). One then sees that 
for the chain located at p* =0, the fixed point at 0* = n 
is always unstable, while that at <jf a = is unstable only 
when k > 4, being stable otherwise. As for the chain at 
p* = 1/a, the point at <f)* Q — is always unstable, while 
that at 4>* = 7r is unstable when k > 4; in both cases, 



de-stabilization occurs via period doubling of the elliptic 
points. One can also make an estimate of the condition 
to be observed for separatrix touching. To do so, let us 
first imagine that we are working in a situation where 
k is small and the dynamics is therefore mostly regular. 
What happens then is that the advance within any par- 
ticular resonant island tends to be slow. One can thus 
approximate p n +i ~ Pn and 4>n+\ — 4>n by their respec- 
tive infinitesimal increments dp and d(f> and finally write 
an expression valid in the vicinity of a m = chain: 
dp/d<j) = — fcsin(0)//(p). Then, with obvious notation, 
one obtains an expression for the separatrix of the lower 
chain: 

t-Pse-p 

/ f{p)dp = k{\+COs{4>sep))- (2) 

The upper separatrix of the (p = 0,4> o = tt) fixed 
point touches the (p Q = l/a,(j> = 0) fixed point of the 
corresponding upper chain (see Fig.(yJ)) when k — k t = 
1/(1 2a 2 ); such value for k is also known as the recon- 
nection threshold in regular regimes. We shall extend 
this definition for the threshold into chaotic regimes. If 
we do that we can draw Fig.(Q). The figure is similar 
to corresponding figures shown in Refs. Jl|Q. It displays 
simultaneously the threshold and period doubling curves 
in the parameter space. The period doubling curve, in 
the m = case analyzed here, is simply a horizontal 
line at k = 4. It is thus seen that for values of a be- 
low a ~ 0.144 period doublings occur before reconnec- 
tion. Previous works have focused interest on the region 
a > 0.144 because in that region one could clearly speak 
in terms of reconnection - recall that in this region re- 
connection takes place before period doublings. In the 
present paper we shall concentrate efforts to see what 
happens deep into chaotic regimes when a < 0.144. 

To start with the investigation, let us consider the 
vicinity of the threshold curve at k = 5. For this value 
of k, elliptic points have been totally destroyed by full 
cascades of period doublings. Therefore, as we increase 
a the theoretical reconnection threshold can be attained 
while the system remains in a deep chaotic regime. Let us 
try to examine how the relevant manifolds of the m = 
resonances behave on the phase-space. The analysis is 
made with help of the panels of Fig. (|^) were we focus at- 
tention on the upper unstable and lower stable manifolds 
(respective orientation indicated by arrows) of the origi- 
nally hyperbolic points (points indicated by black dots) of 
the lower and upper chain respectively; in order to draw 
the manifolds we launch 1000 initial conditions along the 
linearized manifolds, iterating the dynamics forward or 
backwards according to the case. 

First of all we note that in integrable cases separatri- 
ces describe homoclinic loops. Now, even in our noninte- 
grable case, when a is small it is seen that the tendency 
of the unstable manifold of the lower chain is to follow 
the homoclinic loop of the integrable approximation. Of 
course, due to the nonintegrable features the unstable 
manifold eventually starts to execute increasingly large 
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oscillations after it first intersects the stable manifold of 
the same point. As the orbit is about to complete the 
homoclinic loop the oscillations grow and as the oscilla- 
tions grow it may happen that the unstable manifold of 
the lower chain crosses the stable manifold of the upper 
chain as well. But what must be observed here is that 
this latter intersection occurs only after the stable man- 
ifold of the lower chain crosses the unstable manifold of 
the same lower chain many times. In this situation one 
can safely look at the process as closely resembling the 
integrable case, although, as mentioned before, one lies 
in a deeply chaotic regime; this feature is somewhat puz- 
zling and shall be considered in detail later. 

Then, when one increases the value of a, the overall 
topology of manifold crossings appears to undergo a sub- 
stantial change. This can be observed in Fig. (Bp) where 
it is seen that this change in topology is in fact very 
similar to what happens in the purely integrable model; 
to our knowledge this had not been realized before. Here 
the unstable manifold of the hyperbolic point of the lower 
chain makes direct connection with the stable manifold of 
the upper chain. This leads to the opening of a new dif- 
fusive channel connecting the regions located below the 
lower chain and above the upper chain. In other words, 
although one lies in a deep chaotic regime since the orig- 
inal elliptic points have fully period doubled to chaos, 
noticeable changes can be expected as a result of clear 
alterations on the topology of the manifolds. This be- 
havior of the manifolds is rather conservative in the sense 
that while the elliptic points have bifurcated, the man- 
ifolds still try to preserve some aspects of integrability. 
As it appears this feature takes place because near the 
midpoint between the two chains, the orbital frequency 
attains a maximum df(p)/dp = in view of nonmono- 
tonicity. As a result, the local dynamics is relatively lin- 
ear and nonintegrable effects are relatively smaller than 
in other regions of the phase-space. 

We now proceed to show that reconnections in chaotic 
regimes may have direct influence on macroscopic pro- 
cesses, such as diffusion. We actually measure a fraction 
of particles that are transmitted across the region where 
df(p)/dp ~ in a numerical experiment that goes as fol- 
lows. 

A set of 1000 initial conditions is launched in the region 
below the lower reconnecting chain and iterated many 
times. All particles arriving at the region above the upper 
reconnecting chain are reinjected into the initial lower re- 
gion such as to create a steady state in the long run. The 
simulation is optimized by reflecting particles that move 
into the region below the injection momentum p = —2.0. 
The mental picture one can form of the system is that 
of a multicomponent gas placed in a vessel divided by 
a semi-permeable membrane whose role is played by the 
reconnecting chains. The gas is placed below the lower 
chain, and some of the particles are able to crossover up 
to the region above the upper chain. As time evolves, 
one reaches a steady state and we measure transmission 
by computing the number of particles in the upper region 



divided by the total number of particles. 

Figure ([|) shows the fraction of particles transmitted 
after many iterations, versus a, for k = 5 and 6. On 
this range of k, as suggested by Fig.(||), there is no KAM 
curve in the phase space. In Fig.(0) we see an increase in 
particle transmission starting at the reconnection. The 
starting point of the increase can be evaluated analyti- 
cally using Eq.(H). For k = 5 one obtains a = 0.128 in 
good agreement with the simulation. It is also observed 
that for k = 5 the transmission reaches a plateau around 
a = 0.15; presumably at this value the reconnection is 
fully completed. 

In Fig.(^) we finally display surface of sections of 20 
trajectories initially placed at p = —2.0 with uniform 
distribution along <f>. Fig.([|a) is made for a value of a 
prior to the reconnection, a — 0.10, and Fig.(^b) with 
a past the reconnection threshold, a = 0.16. The alter- 
ations in the diffusive pattern suggested by all the previ- 
ous analysis can be seen in those figures as well - while in 
Fig. (^a) the particles remain mostly in the lower region, 
in Fig.^b) particles can be easily transmitted across the 
barrier at p ~ 0. 

To summarize, in this paper we have investigated the 
effect that reconnections involving unstable manifolds of 
deeply chaotic regimes can have on some macroscopically 
observable features like particle diffusion. We have used a 
nonmonotonic map to create reconnecting chains. Exam- 
ining a particular family of fixed point for which m = in 
the notation of the text, we have seen that the topology 
of manifolds of unstable fixed points may have a similar 
behavior as in an integrable approximation. As we vary 
convenient parameters, manifolds of differents chains of 
islands can be clearly seen to undergo reconnection pro- 
cess even in deep chaotic regimes. This is somewhat un- 
expected, since if one lies in chaotic regimes where man- 
ifolds are already strongly interlaced, any reconnection 
would be expected of little influence in macroscopic fea- 
tures. 
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FIGURE LEGENDS 

FIG. 1. Inverse tangent bilurcation and the preceding re- 
connection in a purely integrable case. The figure is con- 
structed with an integrable one-degree-of-freedom Hamilto- 
nian of the type H — p — ap 2 + k cos <j> with (p, (f>) as contin- 
uous canonical variables; k — 5 and a = 0.10 in (a), 0.128 in 
(b), and 0.16 in (c). 

FIG. 2. Parameter space (a, k) for the period one m = 
resonance, and the relevant threshold curves. 

FIG. 3. The reconnection of stable and unstable manifolds 
in a deeply chaotic regime for which k = 5. 

FIG. 4. Transmission fraction as a function of a in the 
deeply chaotic regimes k — 5 and k — 6. 

FIG. 5. Poincare plot of the k = 5 case for a — 0.10 in (a) 
and for a = 0.16 in (b). 
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Fig.1 : Corso and Rizzato 




